Ultracold atoms in [7(2 ) non-Abelian gauge potentials preserving the Landau levels 



O 
(N 



O 



X 



Michele Burrello''^ and Andrea Trombettoni''^ 

'SISSA, Via Bonomea 265, 1-34136, Trieste, Italy 
^INFN, Sezione di Trieste, 1-34127 Trieste, Italy 
(Dated: August 4, 2011) 

We study ultracold atoms subjected to U (2) non-Abelian potentials: we consider gauge potentials having, 
in the Abelian limit, degenerate Landau levels and we then investigate the effect of general homogeneous non- 
Abelian terms. The conditions under which the structure of degenerate Landau levels is preserved are classified 
and discussed. The typical gauge potentials preserving the Landau levels are characterized by a fictitious mag- 
netic field and by an effective spin-orbit interaction, e.g. obtained through the rotation of two-dimensional 
atomic gases coupled with a tripod scheme. The single-particle energy spectrum can be exactly determined for 
bJQi a class of gauge potentials, whose physical implementation is explicitly discussed. The corresponding Landau 

^ ' levels are deformed by the non-Abelian contribution of the potential and their spin degeneracy is split. The 

I related deformed quantum Hall states for fermions and bosons (in the presence of strong intra-species interac- 

tion) are determined far from and at the degeneracy points of the Landau levels. A discussion of the effect of 
the angular momentum is presented, as well as results for U (3) gauge potentials. 

ly^ , PACS numbers: 

■ 

bJO; 

4L ! I. INTRODUCTION 

C ; 

^ ' The recent experimental realization of artificial magnetic and electric fields acting on neutral atoms iH 0] opened the way to 
Q^' simulate ultracold many-body systems in controllable (static) electro-magnetic fields |3]: these synthetic fields can be imple- 
^ mented using spatially dependent optical couplings between internal states of the atoms. This technique has been applied so far 
not only to single-component Bose gases 010, [4], but also to Bose-Einstein condensates with two components isfl^: in isj], 
H using a suitable spatial variation and time dependence of the effective vector potential, a spin-orbit coupling has been realized. 
Besides, the optically-induced effective gauge potential and spin-orbit couplings may be experimentally investigated in the fu- 
ture also in ultracold fermionic gases and Bose-Fermi mixtures. Such a possibility, eventually together with the use of rotating 
traps [f\ and/or the application of more general multipod schemes fs*], envisions the concrete opportunity of manipulating and 
Q ^ studying interacting ultracold systems in a vast class of simulated magnetic and electric fields. 
'— ' i An important perspective motivated by the realization of synthetic magnetic fields is given by the simulation of quantum Hall 
^ ' physics with ultracold atoms. As it is well known, the effect of a rotation on a neutral ultracold gas is equivalent to that of a 
magnetic field on a system of charged particles |9]: since the first experiments with Bose gases in rotating traps 1 10, 11], this 
fact called for the possibility to realize Laughlin states and other quantum Hall states using two-dimensional strongly interacting 
, ultracold gases in rapidly rotating traps (see the reviews ||3l0,[i3.[Iltl)- However the achievement of Laughlin states and quantum 
00 ' Hall regimes with rotating gases turned out to be in general not an easy task, since for typical experimental numbers the rotation 
' frequency ft should be very close to the trap frequency oj 10, [HI]. The recent realization of synthetic magnetic fields using 
I spatially dependent optical couplings adds then new technical possibilities in view of the experimental realization of quantum 
' Hall states. Moreover the first experimental evidences of these strongly-correlated states have been produced for a set of a few 
atoms of ®^Rb loaded in time-modulated optical lattices |T7]. 

Another very active line of research is driven by the possibility to apply controllable artificial magnetic fields on a two- 
(eventually multi-) component ultracold gas, in particular engineering tunable spin-orbit interactions. Spin-orbit coupled Bose- 
Einstein condensates have been experimentally realized fs", o]: in fs*] the spin-orbit coupling had equal Rashba and Dresselhaus 
strengths, but a much wider class of spin-orbit couplings may be realized |8, 15, 16]. At low temperatures, a spin-orbit coupled 
Bose gas can condense into degenerate minima at finite momenta lIT^ : moreover, modifying the interparticle interaction, a spin 
" " stripe phase may appear 1 18]. Equally interesting would be the realization of spin-orbit coupled ultracold fermionic gases: in 
particular, a Fermi superfluid in the presence of a tunable spin-orbit could be studied. The properties of the superconducting state 
in the presence of a spin-orbit interaction lifting the spin degeneracy, with mixed singlet and triplet pairings, include anisotropic 
spin magnetic susceptibility and finite Knight shift at zero temperature 1 19]; such properties have been actively studied in the last 
decade ll20l - i23[l and, very recently, the possible realization and study of their counterpart in atomic Fermi superfluids attracted 
significant attention ll24l - t30ll . 

Spin-orbit couplings acting on a two-component ultracold gas are a particular example of the so-called non-Abelian gauge 
potentials: the single-particle Hamiltonian has in general a term proportional to (p+A)^, where the vector A has non-commuting 
components, e.g. [A^, Ay] ^ 0. It is intended that the vector potential is a matrix (e.g., 2x2 for a two-component gas) and, in 
general, it can have a spatial dependence. In the tripod scheme, three internal quasi-degenerate states are coupled with a fourth 
and the two resulting degenerate dark states are subjected to an effective non-Abelian gauge potential fisi [sTIl . This scheme 
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can be extended to the tetrapod configuration jsTll : a discussion of more general multipod setups is presented in U{2) 
vector potentials acting on ultracold atoms in optical lattices can be implemented using laser assisted tunneling depending on 
the hyperfine levels ll32ll : besides, effective non-Abelian gauge potentials have been also discussed in cavity QED models IBSll . 

Several properties of ultracold atoms in artificial non-Abelian gauge potentials have been recently studied Il34l - l46ll . If a 
general non-Abelian term is added to an Abelian potential exhibiting (eventually degenerate) Landau levels, then the coupling 
between the internal degrees of freedom breaks the degeneracy of the Landau levels: e.g., in lf39ll a matrix generaUzation of 
the Landau gauge has been considered, and the consequent disappearance of the Landau level structure for large non-Abelian 
terms investigated. The possibility of having anomalous quantum Hall effects in suitable artificial gauge potentials has been 
also addressed I4ll - l43ll . The Landau level spectrum in a spatially constant non-Abelian vector potential in planar and spherical 
geometries has been discussed in |45], showing that the adiabatic insertion of a non-Abelian flux in a spin-polarized quantum 
Hall state leads to the formation of charged spin-textures. 

An appealing issue regarding two-dimensional ultracold atomic gases subjected to controllable non-Abelian gauge potentials 
concerns the possibility to use them to simulate and manipulate ground-states having non-Abelian excitations. This interest is 
due to the highly non-trivial topological properties of such correlated states fisl l49ll and to their relevance for the topological 
quantum computation schemes |49J. However, even if the Landau level structure is not broken, in general the excitations of the 
system remain Abelian [42, 44, 46J. In ll46ll . using exact diagonalization, the fractional quantum Hall effect for two-dimensional 
interacting bosons in the presence of a non-Abehan gauge field (in addition to the usual Abelian magnetic field) was studied, 
obtaining that for small non-Abelian fields one has a single internal state quantum Hall system, whereas for stronger fields there 
is a two internal state behaviour (or the complete absence of Hall plateaus). 

In i44ll the Landau levels and the quantum Hall states were studied for a two-component two-dimensional gas subjected to the 
non-Abelian gauge potential = q<7x — -f-y, Ay = qcTy + (where the cr's are the Pauli matrices). The reasons for choosing 
such potential were the following: /) for g = 0, an artificial magnetic field B (perpendicular to the plane xy) is applied to both 
the components and the Landau levels are doubly degenerate; //) for finite q, i.e. for a finite non-Abelian term, the Landau levels 
are preserved, but their degeneracy is lift; Hi) closed analytical expressions for the single-particle energy spectrum can be found 
(see isdH ): iv) last but not least, one can choose the parameters of the laser pulses in a tripod scheme in such way that it can 
be experimentally implementable. For both bosons (in the presence of a strong intra-species interaction) and fermions, explicit 
expressions for the ground-states were obtained, showing that deformed Laughlin states - with Abelian excitations - appear; 
however, ground-states with non-Abelian excitations emerge at the points (e.g., q^ = 3B) in which different Landau levels have 
the same energy. 

The goal of this paper is two-fold. From one side we investigate the effect of general homogeneous non-Abelian terms added 
to the usual Abelian magnetic field (assumed equal for both the components, so that the Landau levels are degenerate in the 
Abelian limit), in order to discuss and classify under which conditions the structure of degenerate Landau levels is preserved. 
From the the other side we present a detailed discussion of the corresponding many-particle quantum Hall states derived in |44l, 
providing additional results and a discussion of the effect of an angular momentum term. We also consider U{3) non-Abelian 
potentials giving rise to lines (in place of points) of degeneracy. 

The plan of the paper is the following: in Section |ll] we study the single-particle Hamiltonian and the Landau levels of a 
two-component two-dimensional Bose gas in an artificial U{2) gauge potential, having doubly degenerate Landau levels in 
the Abelian limit and a non-Abelian SU{2) gauge potential independent on the position. We argue that there are only three 
classes of (quadratic) Hamiltonians preserving the degeneracy of the Landau levels and giving rise to an analytically defined 
Landau level structure, where the eigenstates are expressed as a finite linear combination of the eigenstates of a particle in 
a magnetic field. The first class corresponds to an Abelian U{1) x U{1) gauge potential and it refers to uncoupled internal 
states, whereas the other two are characterized by a truly non-Abelian U{2) gauge potential and correspond to different kinds 
of Jaynes-Cummings models. Then we focus on one of these Jaynes-Cummings classes which is gauge equivalent both to a 
Rashba and to a Dresselhaus spin-orbit coupling: after discussing in Section |III] how this non-Abelian gauge potential can be 
implemented in a rotating tripod scheme, we analyze the Landau level spectrum in SectionHVl where we also present results for 
U{3) gauge potentials acting on a three-component gas. In Section[V]we write the deformed Laughlin states in the presence of 
two-body interactions, while in Section [Vl] we discuss the effect of an angular momentum term. Section [VIll is devoted to the 
study of the ground-states and excitations at the degeneracy points, while our conclusions are in Section [VIII| 

II. SINGLE-PARTICLE HAMILTONIAN AND LANDAU LEVELS 

In this Section we first introduce the single-particle Hamiltonian of a single-component gas in a constant magnetic field in 
order to set the notation for the following results. We then consider a two-component gas characterized by two internal degrees 
of freedom: these two components provide a pseudospin degree of freedom (hereafter denoted simply as spin). We analyze the 
general properties of its spectrum when a non-Abelian SU{2) gauge potential independent on the position, thus characterized 
by a constant Wilson loop |41], is added to a constant magnetic field (equal for both the spins). 

We show that there are only three classes of quadratic Hamiltonians, describing a single particle in such effective U (2) gauge 
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potential, giving rise to an analytically defined Landau level structure (with eigenstates expressed as a finite linear combination 
of the eigenstates of a particle in a magnetic field). The first class corresponds to an Abelian U{1) x U{1) gauge potential 
and it refers to uncoupled spin states, whereas the other two are characterized by a truly non- Abelian U{2) gauge potential and 
correspond to different kinds of Jaynes-Cummings models. 



A. 17(1) gauge potentials 

We first remind the standard case of a spinless atom moving on a plane and subjected to a (fictitious) constant magnetic field 
lIsTll . This case corresponds to a U{1) gauge potential: the Hamiltonian reads 

H = {P, + + {Py + Ayf (1) 

(with units such that h ~ 1 and m = 1/2). The vector potential A in the symmetric gauge is 

B B 

A^^-—y, Ay^—x. (2) 



Introducing the complex coordinate z = x — iy, we define the operators 



and 



^-™(^- + 45.) (3) 



^ (Bz - ids) ■ (4) 



2V2B 

These operators obey the commutation rules 

[D,D^]^1, [L\L]=l, [L,D]^[L,D^=0. (5) 

D and Dt 

can be considered as ladder operators and allow us to express the Hamiltonian (HJ as _ff = 2B [D^ D + 1/2). The 
operators L and enter the definition of the angular momentum of the particle through the relation 

= D^D - W (6) 

such that L and it respectively decrease and increase L^. Each energy eigenstate is degenerate with respect to the angular 
momentum, therefore it is possible to characterize the usual Landau levels by the index n = D. 

In order to maintain the angular momentum degeneracy of the eigenstates, a generic Hamiltonian H for a spinless atom have 
to be independent of the operators L and L^. Imposing this Hamiltonian to be at most quadratic in the momentum we obtain that 
H can be written as a function of the operators D and as 

— = D^D + aD"^ +a*D + bD^'^ + b*D'^ +K,' (7) 

where E' is an overall energy scale, K! is a real parameter, whereas a and h are complex coefficients. It is easy to show that the 
Hamiltonian (|7]i can be recast in the following form 



where the operators V is defined as 



h^et'^t + k: (8) 



(9) 





















a* 1 







in order to satisfy the commutation relation [F, F^] ~ 1, with the complex parameters a and /3 related to the coefficients in (|7]). 
The operators F, F^ allow to express the Hamiltonian in the simple quadratic form (O, which makes evident the Landau level 
structure of the single-particle problem. It is interesting to notice that the mapping from D, to F, F^ corresponds to an affine 
transformation of the space coordinates of the kind: 

i IT If i\ 1/ i\_ As" 
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B. (7(1) X U{1) gauge potentials 

We consider now a system of (non-interacting) atoms characterized by a pseudospin degree of freedom that can assume the 
eigenstates |t) and \l) on the z direction. Therefore, in the above description of the single-particle system, it is necessary to 
introduce also the Pauli matrices a to complete the observable algebra generated by D. D\L and . The Pauli matrices 
and ay couple the two spin components, whereas and the 2x2 identity matrix I describe particles whose states |t) and jj,) 
are decoupled. More generally, to describe a system characterized byaJ7(l) x U{1) symmetry corresponding to two decoupled 
states in the same magnetic field B, the Hamiltonian must be a function of the operators D,D\l and a single linear combination 
of Pauli matrices, k ■ a, that we can relabel as Uz without loss of generality. 

This generic problem is defined by the Hamiltonian 

H = {a^Oz + ao) D'^ + (a^^ + a^) ^ + {^z^z + ^o) ^"^^ + + K) + iK^z + K) D"^ D + XV,. (11) 

Such Hamiltonian describes a particle in the uniform magnetic field B and it is the most general one which is quadratic in the 
momentum, fulfils the U(l)xU(l) gauge symmetry and is independent on the angular momentum. It constitutes a simple gen- 
eralization of the Hamiltonian (|7]i to two non-interacting spin components, thus it can be solved by implementing two different 
coordinate transformations of the kind dTOb for the two inner states. Using the unitary transformation U = Q-i<T^(iaz-ia z)^ j-j^g 
operators D and L)t become 



D^U'^DU = D+J—aaz, ->■ D'U ^ + \ — a*az. (12) 

V i) V 

This unitary transformation, combined with the real space affine transformations JTOb . allows us to define the operator 

r = (ao + aiaz) D + {Po + Picfz) + -fa + -fiaz (13) 

with the parameters chosen in order to satisfy the commutation relation [F, F^] — 1. Apart from constant terms, the Hamiltonian 
(fTTl i can be rewritten as 

Ha = T^T{hol + Kaz)+Maz (14) 

with ho, hz and A4 real parameters. Both the spin dz and the angular momentum Lz are proper quantum numbers and Ai 
constitutes a Zeeman term. Thus, the spectrum of the Hamiltonian is characterized by the presence of Landau levels, degenerate 
with respect to Lz and labelled by both ct^ and n = r^F: if hz and Ai go to zero, all the Landau levels result degenerate also 
with respect to the spin. 

We conclude this Section by observing that every single-particle Hamiltonian of the kind ( fTTT ) can be recast in the form (fT4l ) 
with suitable transformations. Equation (fT4l) makes evident the existence of a quantum number n = F^F which defines the 
Landau level structure. The main characteristic of the Hamiltonian Ha is the fact that it does not depend on a-^ and a-y, and 
therefore it corresponds to a gauge symmetry U{1) x U{1): we will refer to this case as the Abelian limit of the more general 
C/(2) we are going to introduce in the following. 



C. (7(2) gauge potentials 

In this Section we consider particles subjected both to an artificial magnetic field and to a general artificial SU{2) homoge- 
neous non-Abelian gauge potential (simulating a general spin-orbit coupling). Our goal is to characterize the conditions under 
which the non-Abelian gauge potential preserves the Landau levels. More precisely, we want to classify the single-particle 
Hamiltonians with a general homogeneous non-Abelian term added to a magnetic field (equal for both spins) such that the 
following properties are satisfied: 

■pi: Their energy spectrum presents a Landau level structure. 

V2: Every Landau level is degenerate with respect to the angular momentum Lz- 

■P3: In the Abelian limit, the Landau levels become degenerate with respect of the spin degree of freedom. 

The condition WW about the existence of the Landau levels is very general: it is indeed known that for a broad class of spin- 
orbit interactions the spectrum of the single-particle Hamiltonian is composed by eigenstates expressed as infinite series ll52ll . 
For simplicity, in the following, we restrict our attention to Landau levels such that the corresponding wavefunctions can be 
expressed as a finite sum of terms. 
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We observe that the condition l7-'3l is not strictly necessary to obtain a proper Landau level structure, however it is necessary to 
implement the U{2) gauge symmetry which will characterize the Hamiltonian we will investigate in the following. In general, 
terms as the ones in hz and Ai in (fl4l i break this symmetry and do not satisfy the condition ings I Nevertheless they do not spoil 
the main characteristics of the system we will study in section |IV] and can be easily taken into account in the following analysis. 

To obtain the most general Hamiltonians satisfying the conditions 17^11 - IP3I we use the single-particle algebra defined in 
Sections III Al - III Bl Condition l"P2l implies that the spectrum of the Hamiltonian must be degenerate with respect to the angular 
momentum L^- Therefore one obtains [H, L] = [H, Lt] = and the Hamiltonian does not depend on L and Lt, but only 
on D, D\ (t. As in the previous case (fTTT i. spin-orbit couplings can give rise, in general, to terms in the Hamiltonian that are 
not proportional to D'^D. Therefore it is neces sary to introduce a wider class of ladder operators F and F^ generalizing the 
operators D and previously introduced (see p2!] for an accurate description based on the Landau gauge). Condition IP II can 
be rephrased in terms of these generalized ladder operators imposing that there must exist an integral of motion n characterized 
by 

[H,n]=0, [^^,F^]==F^ [n,F] = -F, (15) 

where F obey the commutation relations [F, F^] = 1 and [F, L] — [F, L^] = 0, and there exists an eigenstate 'i'o of F with 
eigenvalue 0: F^Pq = 0. F can be defined as a linear combination of D, D'^ and, eventually, some constant terms. 

The role of the integral of motion n is to label the generalized Landau levels obtained from the couplings between the 
pseudospin states. To satisfy ( fTsT l. n must be chosen in the form 

n = FtF + c-(T, (16) 

where cis a real vector we want to determine. 

We define now the most general single-particle Hamiltonian satisfying the previous conditions, with the constraint that it can 
be at most quadratic in the momentum p (and, therefore, in the operators D and D^). We can divide the Hamiltonian into two 
terms: the first one. Ha, corresponds to the Abelian gauge symmetry U{\) x U{\) and represents the case of uncoupled spin 
components; 

Ha^{ho + Ka,)D^D + M,a,. (17) 
The second term, Hna, is instead the non-Abelian contribution Hna 



Hna = {a ■ a + ao) D'' + {a* ■ a + a^) D+[b - a + boj D'''^+[b* • ct + 6^ j D'^ + {h^a^ + hyGy) D"^ D+M^a^+MyGy 

(18) 

where is a real vector, /i^ is a real vector with spatial part h, and a^, 6^ are complex vectors with spatial parts a and h. The 
total Hamiltonian (apart from constant terms) is given by iJ = Ha + Hna- 
We observe that, if the following conditions are satisfied: 

• a. — e*^"aff, with dji a real vector and Oa a real constant; 

• 6 = e'^^''bj^, with 6^ a real vector and 6i, a real constant; 

• all the vectors A4,h,bji and or are parallel to each other along the direction k; 

then the Hamiltonian can be reduced to the U{1) x U (1) case ( fTTI ) previously studied through proper transformations. In this 
case we have shown that there is a Landau level operator n = F^F commuting with the Hamiltonian. The previous conditions 
imply that the system is not characterized by a proper U{2) gauge potential since only one effective component of the spin, k ■ a, 
enters the Hamiltonian. 

In the generic U{2) case with arbitrary vectors a, b, it is impossible to gauge away all the terms in a • tr. It is therefore 
convenient to search for an integral of motion n defined in terms of the D operators as n — D'' D + c ■ <t; in order to satisfy the 
relation [H, n] = 0, we obtain from the commutation rules the conditions 

e''""/i;c™ = e'™"XiC™ = (19) 
-a„ + 2ie''""a,c^ = , a*„ + 2ie''""a;'c™ = , ao = (20) 
-2bn + 2ie'™"6,c™ = , 26; + 2ie''""&rcm = , 6o = (21) 

Equation (fT9] l requires the vectors h, Ai and c to be parallel; without loss of generality, we can impose them to be in the z 
direction with an appropriate spin rotation and c assumes the form (0, 0, c^). The equations (|20] | and (ISTT i are not compatible in 
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general, unless either a or 6 is zero. Imposing in (l20l l oq = and bo = 0, by using the transformation Equation ( [20b can be 
recast in the form 

3 (a) = 25R (a) x c (22) 

(a) -23 (a) x c. (23) 

The previous equations state that the real and imaginary parts of a have to be orthogonal to each other and orthogonal to c. 
Therefore the condition ( fTSl l implies, in the case 6 = 0, that a must lie in the xy plane and, moreover, that \cz\ — 1/2 in order to 
satisfy ( l22b and ( |23] ). Therefore we can choose a proper spin basis in which 

a = {iq/2,q/2,0) , c= (0,0, 1/2). (24) 

Similarly one can consider the case in which a = 0: to satisfy equation (ISTT i. 6 must lie in the system plane and 3fJ(&) _L 5(&) 
with \cz\ = 1 (which is incompatible with the previous case). 

So far we considered conditions IP 1 l and lP2l and we obtained, in the general U{2) case, that there are two possible Hamiltonian 
classes defined hy b — and a = 0. The transformation (|9]l required to obtain uq — and bo — implies that the Landau level 
operator n has the general form ( fTSl l in terms of the operators F and F^. For the sake of simplicity, hereafter we restrict ourself 
to the case T — D and F^ = 13^, since all the other cases can be studied using the coordinate transformation ( fTOb . Therefore we 
will deal with generalized Landau levels expressed as a finite sum of the eigenstates of D'^D. In this case we are reducing the 
previous Hamiltonians to the following classes: 

• The Jaynes-Cummings class, obtained by imposing 6 = 0, with Hamiltonian 

H ^{E + hzCJz) D^D + MzOz +JC- iqa+D + iqcj-D'^ (25) 

where a± — axiia-y- In the limit hz,A4z — ^ this Hamiltonian satisfies also the condition l7^3l and it is characterized by a 
full U{2) gauge symmetry. Moreover, it can be shown that the Hamiltonian dZSb is gauge equivalent to both a pure Rashba 
spin-orbit coupling and a pure DresseUiaus interaction. This case will be extensively discussed in the next Sections where 
we will show that it can be described in terms of a minimal coupling with a non-Abelian gauge potential. 

• The two-photon Jaynes-Cummings class, obtained by imposing a = 0, corresponding to 

H ^{E + h^(Tz) D^D + MzC^z +IC- iqa+D'^ + iqa^D^'^. (26) 

This Hamiltonian cannot be described by a quadratic minimal coupling with a non-Abelian gauge potential since it 
presents the product between quadratic terms in D and and the a matrices. However it can be exactly solved ifssllsill . 
showing a Landau level structure. Also in this case the condition lP3l is satisfied in the limit hz,Mz — ^ 0. 

Summarizing, up to transformations of the spin basis, there are only three classes of Hamiltonians, quadratic in the momentum, 
that satisfy the conditions VP 1 1 and 17^21 and can be described by generic ladder operators F and F^^. The first one is the Abelian 
class with a U{1) x U{\) gauge symmetry ( fT4l l. The other two are characterized by a full U{2) gauge potential and correspond 
to different Jaynes-Cummings models. In particular we will restrict to the case in which F and F^ correspond to D and ZJt and 
we will focus, in the following Sections, on the class dZSl ) in the limit of U{2) gauge symmetry, fulfilling also the condition l7^3l 



III. ENGINEERING THE NON-ABELIAN GAUGE POTENTIAL 



In this Section we discuss the physical implementation in a rotating tripod system of the non-Abelian gauge potential 

B B 
Ax = qo-x - yyl , Ay = qay + —xl (27) 

(the identity matrix I will be dropped in the following). The corresponding Hamiltonian is 

f B V f B V 

H = Ipx + q(7x - j + \Py + q(7y + I . (28) 

As discussed in the previous Section, the vector potential ( |27] | is representative of a much more general class of single-particle 
systems (1251 ) characterized by the properties WW - fP3\ Moreover the non-Abelian term of dZTl i mimics the effect of a spin-orbit 
couphng and it can be shown to be gauge equivalent both to the Dresselhaus and to the Rashba couphng. 
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The vector potential dZTl ) is constituted by a SU{2) term proportional to the parameter q, quantifying the strength of the 
non-Abelian term, and by the magnetic contribution: A describes a proper U (2) potential whose total effective magnetic field is 



B-2q^ 
B + 2q' 



B = y xA + iAxA=[ ^'^ „ , ) (29) 



where B is proportional to the commutator of the co variant derivatives. It is important to notice that B does not depend on the 
position so that the system is characterized by a translationally invariant Wilson loop as in the cases analyzed in |41]. 

It is well known that the effect of a constant magnetic field can be reproduced in a rotating frame thanks to the Coriolis force 
(see, for example, |7]); at variance, the SU{2) contribution can be obtained through proper optical couplings in a system of 
atoms showing quasi-degenerate ground-states as described in ifTsi Issll . Therefore, to engineer the effective gauge potential 
A{q) (l27T i. we will consider a rotating system of the so-called tripod atoms whose coupling is described in llisll and depends on 
the Rabi frequencies VLi 

Oi = i7sin(6')cos((/))e'^i , = sin (6*) sin (0) rjg = 12 cos (6*) e*^^ (30) 

where Si and 5*2 are functions of the position and of the parameter q, while the angles and 6 and 6*3 are chosen constants. These 
frequencies describe the couplings of three quasi-degenerate ground-states, characterized by different hyperfine levels, with an 
excited state. This interaction give rise to two different dark states whose dynamics is described by the effective Hamiltonian 
(l28l l with vector potential given by ( |27] |. Such dark states constitute the different pseudospin component |i) coupled by the 
non-Abelian term in H. 

In the following we determine the dependence of the frequencies Hi on q to obtain, through a proper gauge transformation, 
the Hamiltonian (l28T l for the atomic gas in a rotating frame. 

Let us consider first a system of tripod atoms in an inertial frame of reference characterized by a non-Abelian gauge potential 

Ji, a scalar potential Kot = '^[A.) + as the one described in ifTsIl and a harmonic confining potential ujr^ lA where = x^+y^ 
(notice that we are working in units in which m = 1/2). The corresponding Hamiltonian reads 

i/iF = (p + Xj' + ic'V^+Kot. (31) 
Once the whole system is put in rotation with angular velocity il, the Hamiltonian in the rotating frame of reference reads 13; 

H,ot - (p + X) ' + ^2 + nU + ^rot , (32) 



where we introduced the gauge invariant angular momentum 



L = r X 



(p + 1) (33) 



and all the coordinates are now considered in the rotating frame. It is useful to rewrite iJ^ot introducing the gauge potential 



B B 
Ax = Ax--y, Ay^ Ay + -X (34) 



where we imposed B = uj. We obtain 

H,ot={p + A] - AL, + Kot (35) 



with A = UJ — fl. 

Our aim is to identify the correct family of Rabi frequencies, Vmt and gauge transformations such that the Hamiltonian Hrot 
can be cast in the form 

Hl= (p + aY -AL, (36) 

with A given by ( |27] l and — r y. p being the usual angular momentum in the rotating frame. As we will show in section FVTl 
iJi can be exactly solved and, in the hmit — > w, it becomes the Hamiltonian (|28l l. In particular we need to have 

A = {qax,q<Ty) , (37) 
Vrot = qA{xay - ya^)- (38) 
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In order to obtain, in the rotating frame, the potential A in dZTl l starting from the potential A 

An = - cos^ ,/)VS'23 - sin^ (/)VS'i3 (39) 

A22 = - COS^ e (COS^ ^V^ig + Sin^ (/)VS'23) (40) 

A12 = -cos^ Qsin20VS'i2 -iV0^ (41) 

given in fl3\, we need a suitable unitary gauge transformation 0(r) f53\. In particular the field transforms as ^ — > OAO^ — 
iOVO^ and thus we must have 

OAO^ - iOVO"^ = 1= [qa^ , qay) . (42) 

From the definition of A one can see that, choosing a constant (f), it is not possible to obtain Ay cx (jy, but it is possible to 
check that we can obtain Ay = kl — qaz and Ax — q<Jx for a suitable choice of the parameters as functions of q. Therefore the 
gauge transformation we apply is 

^ , with O = e^'^'^-^f (43) 

with k to be defined in the following. In this way we obtain 

■Ax = q(Tx ^ Ax^ qax , Ay = kl - qa^ Ay = qay. (44) 



We have also to consider that the scalar potential in lITsIl is affected by O as O + <!>) O^: then, in order to obtain (|36] | out of 
, we need to have 



O + $) Ot ^ Kot = Af X i (45) 

and then 

V + <^> = -qA{ycrx+xaz). (46) 

We are now in position to find the suitable parameters satisfying Equations (|44] | and (|46] |. First of all we impose = 7r/4 and 
53 — cost; then, from the definition of A we obtain that 

dx (5i + ^2) = (47) 
-cosedxiSi- S2) = 2q (48) 
2k-2q = -dy {Si + S2) (49) 
2g + 2fc = - cos^ edy {Si + S2) ■ (50) 

A possible solution is given by 

Si^X{x + y), S2^X{-x + y) (51) 
with A = —q/ cos6. From the last two equations we obtain 

cos^ 6* - 2 cos 6* - 1 = ^ cos 61 = 1 - V2. (52) 

It follows that 

q q 1 + COS^ 9 r- 

A = ^ = -J^, q = -^q. (53) 

COS t' v2 - 1 2 cos 

The corresponding Rabi frequencies are 



^i = n'e"'^^ , ^2 = Q'e"^^^ , ^3 = -i7'V V2 - le'^^ ~ -OMn'e''^^ (54) 
with i?' and 5*3 arbitrary, so that the right A are obtained after the gauge transformation. 
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Let us consider now the scalar potentials; imposing = 7r/4 and cos 6 = 1 — \pl, we find from ifTsIl and from ( |46] | 

+ = .!l±ii + sin^ 61 = -gAx (55) 

V'22 + $22 = cos^ e + V3 sin^ 6* + cos^ d sin^ 6I = qAx (56) 

V\ - V2 

Vi2 + $i2 = ^ cos9 = -qAy. (57) 

The solution is given by 

Vi = -qAx - XAy - \^ sin^ e (58) 
V2 = -qAx + XAy - siri^ e (59) 
V3 = V2qAx (60) 

with A given by ( 1531 ). 

This choice for the scalar potentials Vi completes the set of parameters ( |54T ) needed to obtain the Hamiltonian ( |36] ) in the 
rotating frame. We notice that it is possible to modify the previous derivation of the scalar potential in order to obtain a Zeeman 
splitting reproducing the term proportional to A^^ in ( |25] l. 

We conclude this Section observing that it is impossible to obtain the gauge potential A defined in Equation ( |27T i (or a gauge 
equivalent version) using only the gauge potential A (I39I40I41| | defined in 1 15] without introducing other physical elements such 
as the rotation of the system. Indeed, applying the gauge transformation i43[ we can express A as 



B . B 

"2 



Ax = -—y + qax, Ay = —x - qa^ + k. (61) 



This form of the potential A is real and does not depend on ay. Therefore, imposing A = A and considering Equation ( |4TI ). we 
obtain that the parameter (p entering the definition of the Rabi frequencies ( l30b must be independent on the position, otherwise 
an imaginary term proportional to V0 would appear Let us consider now the term An: from the equation ( [39] l one obtains 
that dyAn^x — dxAn.y. However this is not the case for An^x and An^y in ( |6TI ) unless B = Q. Therefore to obtain the term 
proportional to the magnetic field B one needs either more complicated multipod schemes or an additional physical mechanism 
(in our analysis the rotation) besides the construction of the non-Abelian gauge potentials for tripod atoms. 

IV. LANDAU LEVEL SPECTRUM 

In this Section we discuss the diagonalization of the single-particle Hamiltonian in the non-Abelian gauge potential described 
in the previous Section: for the sake of simplicity, we will begin our analysis studying the Hamiltonian ( |28] ). corresponding to 
the limit — > cj, which is necessary to satisfy the condition l7^2l in section HTCl In the next Section, we will address the case in 
which also a term linear in the angular momentum ( |36] | is present. Two-body interactions will be introduced in Sections IVl and 
[VAl 

As discussed in Section|IIl the Hamiltonian ( |28] ) can be decomposed into two terms: the Abelian one. Ha, and the non-Abelian 
one, Hna ■ One has 

Ha = 2q^ + B + ^{Bz- Ad,) {Bz + 4d,) =2q^ + B + ^dU (62) 
n ( ^ -iBz-AidA fO -id\ 

where we defined the operator 

d = Bz + 4,di = 2V2BD (64) 

proportional to the previously defined operator D so that the relations [d, z] = and [d, d^] = 8B hold. 

Introducing the standard gaussian wavefunction = e~'^ , one has d^o — and d^ipo = 2Bzipo', thus, for Ha, we obtain 
the usual Landau level structure of the eigenstates pill , degenerate with respect to the angular momentum n — m: 

V'n,™ = X A- « D^'-L-^i'o- (65) 
(SB)- ^ 
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The corresponding energy levels are 



E,, = 2g2 + 2B ( n + i ) . 



(66) 



The eigenvalues of Hna are — ±2qy/2Bn and the corresponding eigenstates (p^m can be expressed in terms of the eigen- 
states of Hn as 



"071— l,m It) ± i'n^m It) 



where the following relation holds for n > 1: 



2V2Bnilj„ 



(67) 



(68) 



In general the non-Abelian term mixes the {n — 1)"^ and the n*'' Landau levels; but there are also uncoupled eigenstates 
fQ = "00. m It) with eigenvalue A = for every -00. m in the lowest Landau level. The spectrum of Hna is similar to the 
one obtained in the relativistic case typical of the graphene systems i56ll : in particular, these results are analogous to the ones 
obtained in |42 ] starting from the Dirac equation in an anisotropic regime, and we can notice that Hna corresponds to the known 
Jaynes — Cummings model, as discussed in SectionHJ 

We can now diagonahze the whole Hamiltonian using as a basis the functions ipn ■ it is 



Hip^ = i2q^ + 2Bn ± 2qV2B^j ifit ~ Bcp^, 
so that, for n > 1, the Hamiltonian is splitted in blocks Hn of the form 



For the uncoupled states one has 



The eigenvalues of H are therefore 



2g2 + 2Bji + 2qV2Bn -B 

-B 2q^ + 2Bn - 2qV2Bn 



H^o = (B + 2q^)^o- 



et = 2Bn + 2q^ ± ./WTWBn 

and its (unnormalized) eigenstates are 

Xtr,^ =[b + 2qV2B^ T VB^ + Sq^Br^j ^n-i.n It) + {b - 2qV2B^ ± ^B^ + Sq^Bn^ V„,„. |t) , 
where we made explicit the angular momentum degeneracy. A plot of the eigenvalues of H is presented in Fig. [T| 



(69) 



(70) 



(71) 



(72) 



(73) 




q /B 



FIG. 1: Energies of Xn+i (black solid line) and xt (red dashed line) for n = 0, . . . , 10 as a function of q^/B. The crossings of the 
ground-states occur in the degeneracy points g^/B = (1 + 2n), denoted by solid circles. 
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We notice that the angular momentum is not a good quantum number for this states unless we introduce also a spin-1/2 
component defining a total angular momentum J — L + S. J commutes with both the Hamiltonians ( |27] l and ( [34] l and 
JXnm = [n — m — 1/2) Xn,m- Moreover, the wavefunctions xt,m ^re eigenstates for the operator n = D + (7^/2 defined 
in ( fT6] i: thus, they constitutes deformed Landau levels. 

Analyzing the spectrum, we can notice that there is a correspondence between the usual Landau levels and the states x's ll44ll : 
every Landau level is splitted into two parts corresponding to the states Xn-i ^^'^ Xn ^i^^, in the case g — > 0, their energy 
becomes approximately E„~i ± 4:q^n. Therefore, for q — > one recovers the usual Landau levels structure characterized by the 
(double) spin degeneracy as prescribed by the condition lP3l The non-Abelian term of the Hamiltonian removes this degeneracy 
through the coupling between the {n — 1)*'' level with spin up and the n}^ with spin down. The deformed Landau levels for 
q > can be defined also considering the Landau gauge; in 114611 it is shown that, in this case, the eigenstates are distinguished 
by the Z2 symmetry obtained by the parity transformation x — >■ —x. 

Varying the value of the parameter q^/B, measuring the ratio of the Abelian and non-Abehan contribution in the gauge 
potential, the eigenvalues show an interesting pattern of crossing points (see Fig. |2|i: each pair of eigenstates of the kind Xa 
and Xb becomes degenerate for 

^ = ^ (a + 6 + VI + 4a6) (74) 

and the energy of the crossing is Ec (a, b) = {a + b) B. Instead, a pair of different eigenstates of the kind xt ™d Xb has a 
crossing only if a < fe; in this case the degeneracy point is 

^ = i (a + 5- VTTI^) : (75) 

also for these levels the corresponding energy is £c (a, 6) = (a + b) B. Therefore all the energy level crossings are characterized 
by an integer energy in units of B. 



SSI '° ° ° ° 



FIG. 2: Crossing points of the eigenvalues: black (gray) circles indicate the crossings between states Xn+i corresponding to the ground-state 
(excited states). Whereas circles denote the crossings of the states Xn+i the states xt, (n = 0, . . . , 10). The coordinate / B for these 
degeneracy points are provided in Equations ( I74I75I ). 

As shown in Fig. [T]the uncoupled state family (fTTl i. corresponding to Xa^ is characterized by the energy = B + Iq^ , which 
is higher than the energy ^2B + 2g2 - y^WTSB^ of Xi- Therefore Xi is the ground-state family of the system for 
q^ < 3B (the general case with q^ > 3B is analyzed in section [VA] l. We can rewrite each state of the family Xi in the form 

xr = e-Tl-l^ (c^ ,p II) + 2c^,iBzP \i) - Aci,AP \i)) , (76) 

where P is a generic polynomial in z and we defined the constants 

Cf,„ = B + 2qV2B^ + ^ B"^ + ^q^Bn (77) 

c^,„ = i(^B- 2qV2B^ - ^/ B^ + Sq^Br^ {2^23^^ ^ . (78) 

It is also convenient to introduce the operator 

Qi = c^.i<7^ + ci^id^ (79) 
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that allows us to map uncoupled states in xt states in Xi Thus, we can rewrite ( l76b using the operator Qi 

Xr -^1 (P(z)e-*l-l'|i)) . (80) 

Due to the degeneracy of these ground-states, one can build also wavefunctions minimizing supplementary terms in the 
Hamiltonian (l28b : for instance, we can introduce a repulsive potential for the spin up component of the form (z, Q = 
5{z — C,) It) (tl where C, plays the role of the coordinates of a quasi-hole in the spin up wavefunction. The corresponding 
single-particle ground-states are 

<j,{zX) = Qi({z~C)Q{z)e-^\^\"\\)) (81) 

where Q{z) is a generic polynomial. We observe that the wavefuntion density for the spin up component in (f) goes to zero for 
z = C, while the spin down density, in general, does not, due to the derivative term in We can also consider a repulsive 
potential not affected by the spin, V {C,) = 5 {z — Q: in this case we obtain as a ground-state the wavefunctions 

$ (z, C) - Gi [{z - Cf Q{z)e-^\'\" i;)) (82) 

having a vanishing density in ^ both for the spin up and the spin down component. Notice that it is impossible to create inside 
the space Xi ^ wavefunction with a zero spin-down density in C, and a non- vanishing spin up component. 

With respect to the Hamiltonian (|27] |. these excitations are gapless; however they increase the total angular momentum of the 
system, and, as we will show in section [VTl this implies an increment in energy once we consider the case A^w — r^T^Oin 
dMll. 



A. U (3) non-Abelian gauge potentials 

As shown in iIsHIstIi . it is possible to engineer gauge potentials involving a higher number of internal states. For instance, 
considering atoms with a tetrapod electronic structure, one can obtain three degenerate dark states, which we denote by |+), 
|0) and |— ): this corresponds to an effective spin 1 and it allows to mimic the effect of an external C/(3) non-Abelian gauge 
potential. 

In particular we can generalize the construction of the previous Section to the following potential: 

^ /O a 0\ ^ / Q -ia Q \ 

A, = a /? , Ay^-xl+\ia -z/? (83) 

^ \0/30y ^ \{) ip Q j 

where I is the 3x3 identity matrix, B is the effective magnetic field and a and /3 are two arbitrary parameters giving the 
coefficients of different Gell-Mann matrices. Ax and Ay do not commute and they describe a particular family of effective 
homogeneous non-Abelian C/(3) potentials whose corresponding non-Abelian magnetic field is 

/B-2a^ \ 

B=\ S + 2a2-2^2 (84) 
\ 3 + 2(3"^) 

which is translationally invariant as in Equation (|29] l. These potentials are similar to the one chosen in iIstIi to simulate Weyl 
fermions through multi-component ultracold atoms in optical lattices. 

Given the potential ( [83] l, the minimal coupling Hamiltonian assumes the form: 

/ 2a^ -iad \ 
H{a,(3) = B + -dU+\iaS 2a^ + 2p^ -i^d] (85) 

4 V i0 2/32 J 

where we used the operators d and d^ defined in (|64] |. The first term in H (a, /3) is the Abelian term proportional to the identity, 
while the second one describes the non-Abelian interaction, depending on the parameters a and /3, coupling subsequent Landau 
levels with different spin as in the U{2) case. Therefore, for each n > 2, we can identify three families of eigenstates obtained 
by hnear superpositions of the states ipn,m |— ), ipn-i,m |0) and ^n-2,m \+), where m is related to the angular momentum. In 
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particular, given n > 2, the corresponding eigenenergies of H (a, /?) are the solutions e„ of the following eigenvalues equation 
(see Fig. [3]for the case a — (3): 

/2a^ + B{2n-3) -i2a^2B [n-X] \ 

i2ay/2B {n - 1) 2a^ + 2/3'^ + B (2n - 1) -^2/3^/2B^ e„I. (86) 

V i2pV2B^ 213^ + B{2n+l)J 

Like the case of the Jaynes-Cummings coupling, the spin degeneracy of the Landau levels is removed and the eigenstates of 
(ISST i are also eigenstates of the total angular momentum J. 

In analogy with the previously discussed U{2) potential, there are also other eigenstates of the Hamiltonian dSSl l corresponding 
to the uncoupled states V'o,m |— ) with energy Eq = 20^ + B and to the family of (unnormalized) "doublet states" defined by 




with energy 

ef = 2B + + ^2 ± ^J{a^ - Bf + 8BP^. (88) 

We observe that in the limit a — > or /3 — > the results of the previous Section are recovered. In fact, if either a or /3 goes 
to zero, the resulting gauge potential is an effective U (2) potential of the kind dZTl l and one of the spin states remains decoupled 
with respect to the others. 




11 > \ > \ > \ I \ > \ > I 

0.5 1 1.5 2 2.5 3 

a = p 



FIG. 3: The energy levels of the Hamiltonian l l85t are plotted imposing the constraint a = (3 (with B — 1). The black lines represent the 
eigenvalues En iS6\ . the purple dashed lines represent the energy of the doublet states and the red line represents the energy eo of the 
uncoupled states. One can see that, in the limit a = /3 — )■ 0, the Landau level energies are recovered with a threefold pseudospin degeneracy. 

Finally Fig. |4]shows that it is possible to recover a triple degeneracy of the ground-state for particular values of the parameters 

a and /?. The figure shows a triple degeneracy occurring for B — 1, a ^ ^^(9 + \/73) /6 and (3 — •\/2/3 between the uncoupled 

state and two eigenvectors of (l86T l obtained for n = 2 and n = 3 at the energy e ~ ^B. We also checked that lines of doubly 
degenerate ground-states occur in the plane defined by a, (3 (single points with triple degeneracy belong of course to these lines). 
Further details on the spectrum of Hamiltonian dSSl ) will be presented elsewhere. 

V. TWO-BODY INTERACTIONS AND DEFORMED LAUGHLIN STATES 

In the previous Section we described a single particle with two internal degrees of freedom in the non-Abelian potential 
(IZTT i: we consider in this Section a system of N atoms, introducing two-body repulsive interactions. Denoting by gi the 
(dimensionless) scattering length between particles in the same internal state and by go the scattering length between particles 
in different internal states, we can write the interaction Hamiltonian as 

N 

Hi^J2{gin,+gono)5{z,-Zj). (89) 
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FIG. 4; The first energy levels of the Hamiltonian l l85t are plotted as a function of a for /? — and B = 1. The black lines represent the 

eigenvalues e„ ( I86t , the purple dashed line represents the energy of the first doublet states and the red line represents the energy eo = 7/3 

of the uncoupled states. For this choice of the parameters a triple degeneracy of the ground states appear for a — \ {% -\- \/73) /6 — 1.71. 



Here Hi is the projector over the space in which the particles i and j have parallel spin states (|tt) or lii))> whereas Ho is the 
projector over the space in which i and j have antiparallel spins (|ti) or |it))- We will consider both bosonic and fermionic 
gases, keeping in mind that for fermions it is 51 = iflolfTTll . 

An arbitrary two-particle state in which both atoms are in Xi can be described as 



* = ^1,1^1,2^ (21, Z2) e-s(l-il^+l-^l^)/4 1;;) 



(90) 



where Qi^i, defined in (|79] |. refers to the coordinate Zi, and P is generic polynomial in zi and Z2- With vanishing inter- 
species interaction {g^ — 0) and strong intra-species interaction, ^ has a zero interaction energy if its components |tt) > lii) 
vanish when zi — > Z2'- for fermions, this is assured by the Pauli principle; whereas for bosons the strong intra-species regime 
corresponds to gi ^ B,q and the two-body wavefunction ^ has to fulfil the requirements 



(91) 



Every antisymmetric polynomial P (zi, Z2) = —P {z2, zi) obviously satisfies these constraints, and, in general, all the fermionic 
functions ^ {zi, Z2) guarantee that the intra-species interaction gives a zero energy contribution. 

If we add also an inter-species repulsive interaction, such that go B,q, the two-particle wavefunction (|90] | must satisfy the 
further constraints 



= 



(92) 



in order to be a ground-state of Hj. These relations hold, for instance, in the case P — {zi — Z2)™ with m > 1. In the case 
m = 2 the inter-species interaction is zero, but not the intra-species one, whereas for m > 3 every repulsive potential Hj gives 
a null contribution. In the following we consider the regime given by go = and (for bosons) gi ^ B > 3q^. Under these 
conditions we can generalize the previous results for the case of N atoms. 

The fermionic states, antisymmetric by the exchange of every pair of atoms, have a zero interaction energy; thus, a possible 



ground-state of the A^-particle Hamiltonian H = J2k 
single-particle energy, is given by 



N 



Hk + Hj, with all the atoms in the Xi space in order to minimize the 



N 



A [Po (^1) , ■•■ , Pn-1 [zn)] e 



(93) 



where A is the antisymmetrization over all the coordinates Zi and P,„ are different polynomials. Generalizing this kind of many- 
body states, it is easy to define a deformation, due to the non-Abelian potential, of the common Laughlin states. If we choose 
Pm{z) — 2™ with m — 0, — 1, we obtain the usual Jastrow factor A [Po...Pn-i] = Y[i<j i^i ~ ^j)- More in general, 
given a Laughlin wavefunction 



N 



A 



(m) 
N 



lu - ;) 



(94) 
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with TO odd, the state 

N 

^^"'^ ^YlGi.jA^N^ (95) 

j 

is a ground-state of the Hamiltonian H: every atom lies in a superposition of states Xi and the antisymmetric wavefunction 
causes the intra-species interaction energy to be zero. 

Also for bosons, symmetric under the exchange of two particles, there are states that have a zero intra-species interaction. For 
instance we can consider for an even value of to > 4. In this case for each pair of particles with — )■ Zj the wavefunction 
vanishes at least as (z; — Zj)'^, thus the interaction energy (and also its inter-species contribution if go 7^ 0) is vanishing. 

Therefore it is important to notice that the introduction of the non-Abelian gauge potential (in the regime < 3B) implies 
that the highest density deformed Laughlin state with null interaction energy has a filling factor 1 /4 instead of the usual filling 
factor 1/2 that characterizes systems of rotating bosons ItI IssIi with a contact interaction. Thus we expect that the introduction 
of the SU{2) potential gives rise to the incompressible state as numerically observed for small values of the chemical 

potential in the weak-interacting regime |46]. Such state is absent in the case of a pure magnetic field and it can be considered 
as a signature of the effect of the potential ( |27] i. 

The state describes in general an incompressible fluid of spin-1/2 particles, as it can be shown calculating its norm; one 
finds 



N 



a(™) 



(|c,.ir + 8i3|c,,i|^)''(A(r^|A(r'), (96) 



where we considered that all the single-particle states involved in are in the lowest Landau level and therefore are eigen- 
states of dS . Thus the norm / can be easily written in terms of the one of the Laughlin state, and one can apply the argument 
also to ^-t"). 

This is true also if we consider quasi-holes in the Laughlin state as, for example 



N / N 



=n^M (^n(^.-ci)'(^.-c2)'j a(,"\ (97) 

since each atom in the Laughlin state is in the lowest Landau level, and thus the operators Qi modify only the norm of the states 
by a constant factor for each atom. 

This correspondence highlights the nature of these excitations since it allows us to state that the Berry phase due to the 
adiabatic exchange of the pair of quasi-holes Qi and ^2 is the same of the one characterizing the corresponding quasi-holes in 
a simple Laughlin state. This kind of excitations are therefore Abelian anyons and their braiding statistics is ruled by the usual 
exchange properties of the Abelian states in the fractional quantum Hall effect Issll . 



A. Generalization to higher value of /B 

So far we referred to the case < 3B in which the ground-state family is provided by wavefunctions of the kind Xi- 
However, as shown in Fig. [T] for higher values of q different ground-states families alternate. Therefore it is necessary to 
generalize the previous results also for > 3B by defining the family of operators Qn describing the deformed Landau levels 
of the kind x - From ( |74| ) one sees that Xn>i '^^e ground-state family for 

{2n-l)B <q^ <{2n+l)B: (98) 

its energy varies, in this range of q, from {2n — 1) B to {2n + 1) B. To describe the ground-state family Xn we generalize 
the operator Qi ( |79l ) introducing the operators 

where the constants „ and c^_„ are defined in ( I77I78I ). The ground-state wavefunctions can be expressed in the form 

X- = Gn (P(z)e-^l^l' 1;)) . (100) 
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Using these expressions and following the procedure shown in the case of it is possible to obtain the appropriate many- 
body wavefunctions for each value of B and q, with n chosen in order to satisfy ( |98l ) (the case of the degeneracy points q^ = 
{2n + 1) B will be analyzed in section [VTlb . In particular, for an arbitrary n, all the antisymmetric states given by 

N 

*i'"^=n^"-j^ir\ (101) 

with A^^-* an odd Laughlin state ( |94l ). are fermionic ground-states unaffected by repulsive intra-species contact interactions 
(here G,i,j indicates the operator applied to the atom j). 

For bosons having repulsive delta interactions (both intra-species and inter-species) one has to consider the derivatives present 
in the operators Qn to find a ground-state having zero interaction energy. The highest order derivative in Qn is given by the 
term d" in the \l) component of each particle. Therefore, to identify the smallest even power m in ( llOlb annihilating a delta 
interaction, one has to consider for each pair of particles the component: in order to make the repulsive delta interaction 
null there must be a factor (zj — Zj)™ with m > 2n. Once the polynomial order of the wavefunction (IIOII ) is high enough to 
make vanish the interaction in the HI) component, then also all the other components give a null contribution. Therefore the 
bosonic wavefunction ^1'^^"'''^' is the ground-state with the smallest polynomial order for generic repulsive delta interaction. 

As a consequence, ^P^^""''^'' defines the maximum filling factor i^n — ^/ {2n + 2) over which the interaction energy among 
bosons cannot be zero. In the range ( |98] l only states with v < 1/ {2n + 2) can have a null interaction energy for Zi — ?> zj in 
the component. This result is consistent with the numerical data obtained in ll46ll where it is shown that, in the case of 

q^ < 3B, the Laughlin state with 1^ — 1/2 has a positive interaction energy, whereas for = 1/4 the energy is exactly zero. 
In this regime, in fact, m must be at least 4 to give a true ground-state, as already observed in the previous Section. Moreover, 
the maximum filling factor decreases as q^ / B increases, therefore the role of interactions becomes more and more important 
if we consider higher filling factors for high values of q. This could explain why, for i/ > 1/2 > Vn, there are no numerical 
evidences of incompressible bosonic states after the first Landau level crossing (q^ > 3B) ll46ll . 



VI. THE EFFECT OF ANGULAR MOMENTUM 



As mentioned in Section |III1 the single-particle Hamiltonian ( |28] l can be considered as the limit of the Hamiltonian iJ^ 
when the angular velocity ft approaches the trapping frequency w. However, so far we considered only cases in which the 
condition 17^2 1 in Section HTCl holds and we neglected an eventual energy contribution of the angular momentum. In this Section 
we analyze the effect of the angular momentum term in the Hamiltonian ( |36] l in order to calculate the energy of the states 
(and their excitations) and to determine the constraints which A = u — fl must satisfy not to spoil the Landau level description. 

Let us consider the single-particle Hamiltonian ( 



Hr 



P + A 



The term proportional to A is spin-independent and it does not affect the non-Abelian contribution Hna in ( l63] l. Using its 
eigenstates „ (67\ for n > 1 we can split Hl into blocks of the form 



HL,n,m = "^q"^ + 2Bn - A ( n - to - i 



The eigenenergies of are therefore 



f2qV2B^ -B+% 

A 
2 



V-B + A -2qV2B^ 



2Bn + 2q^ - A{ n — m 



± 



B- 



Sq'^Bn 



(102) 



and the corresponding (unnormalized) eigenstates are 




lt> + 



A 



B - y - 2qV2Bn ± 




Ii>- (103) 
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We see that the uncoupled eigenstates are unaffected by the angular momentum term. 

Since the coefficients in the definition of the eigenstates do not depend on m, then one can redefine the constants c-|- „, cj^ „ 
and the operators Qn independently on m. Moreover, the Landau level structure holds also for Hl (condition IhD, and the 
Landau levels are energetically distinguishable provided that the energy contribution of the angular momentum remains small 
with respect to the Landau level spacing. In fact, all the energy levels (1102b have a term which is linear in the total angular 
momentum J = n — to— 1/2 (increasing their energy if J < 0). Thus the states with lower values of to are favoured: this 
implies an energy gap A for the creation of single-particle quasi- holes of the kind (HB and 2 A for ^ (for C = 0). 

In order to understand the stability of the Landau level structure, and thus of the deformed Laughlin states, we have to consider 
a multiparticle wavefunction describing N atoms in a ground-state of the kind (llOlb for 2n — 1 < / B < 2n + 1. In 
this case a particle corresponding to the highest value of the modulus of the angular momentum term \J\ acquires an additional 
energy ^ NmA that must be smaller than the gap with the next deformed Landau level. This gap, for values of far from the 
degeneracy points, can be calculated evaluating the energy difference of x~ and Xn-i ™ '^he crossing point between Xn-i and 
Xn+i ^ii'^ i'^ be approximated by B/{2n). Therefore the Landau level description of the multiparticle states remains accurate 
if A < B/{2nNm). 

Let us analyze more in detail the regime characterized by small values of A. The deformed Laughlin states ^1™^ have to be 



defined with the appropriate corrections in the operators Qn ( l99] l since the constants c-\ and cj. must include A coherently with 
(1103b . The corresponding energy contribution of the total angular momentum turns to be for ^l™'' 

_ AJ*,(;™) = AiV _ m - n + « l^rn^^^^\ (104) 

Therefore, the Laughlin states with smaller m (and higher density) are energetically favoured: thus, in the case of fermions, far 
from the degeneracy points, the ground-states are described by the v — 1 filling factor states 5*^, whereas in the case of bosons 
the ground-states are of the form VE'^""''^. 
A quasi-hole in the position C = 

N N 

J i 

acquires an additional energy iVAfc with respect to the corresponding ground-state, because it changes the angular momentum 
of each particle by k. This energy can be considered the gap for the creation of a quasi-hole. 



VII. DEGENERACY POINTS AND NON-ABELIAN ANYONS 



An important property of the single-particle spectrum of the Hamiltonian ( |28] l is the existence of degeneracy points corre- 
sponding to the values q^ ~ (1 + 2n) B: this makes possible the occurrence of many-particle ground-states having non-Abelian 
excitations lliill . In these degeneracy points the two lowest energy levels cross (possibly generating a first-order phase transi- 
tion 1 46]) and the ground-state degeneracy of the single-particle is doubled. In these points an atom in the (non-interacting) 
ground-state can be described by all the superpositions of wavefunctions in Xn ^i^^ Xn+i- However, if we consider the angular 
momentum term in the Hamiltonian (see the previous Section), then the states with a lower angular energy are favoured and the 
variation of the parameter q around the degeneracy points gives rise to a crossover, as we will describe in lVllCI 

When the intra-species interaction between atoms is introduced, the doubled degeneracy of the single-particle states implies 
a novel form of the multiparticle ground-state which is quite different from (193195b . We first consider the case of the first 
degeneracy point, q^ = 3B, in which there is the crossing between a ground-state of particles in Xi and X2 (o^r conclusions 
will be later extended to all the other degeneracy points). For the single particle the basis of ground-states is defined by the 
set {Giz'^ipo (z) , G2z"^ipo {z) , with m G N} and, analogously to the previous sections, we have to distinguish the case of 
interacting bosons and the one of (free) fermions. 

A. Fermionic gases 

We first analyze the fermionic case: the highest density ground-state function of the Hamiltonian H for 2N atoms is given by 

^c^A [Gi^bo, Giz^Po, ... , ^iz^- Vo, ^2^0, 52z4>o, ... , ^2^"^" Vo] \U ... ;> , (106) 
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where A implements the full antisymmetrization over all the atoms. Because of the double degeneracy, the wavefunction ilc 
describes an atomic gas with filling factor ly — 2. The state ftc is obtained through the Slater determinant of the single-particle 
wavefunctions with the lowest angular momenta | J|, up to the power z^"^. Therefore, considering also the angular momentum 
contribution in the Hamiltonian, it is the true ground-state for the system. 

The double degeneracy makes ilc very different from the case ( |93] i: in the degeneracy points the antisymmetrization hides 
a clustering of the particles into two sets of N atoms, say A and B, that are physically different and refer to states in Xi and 
in X2 ■ The two clusters must have the same number of atoms in order to minimize the contribution to the energy given by the 
angular momentum. The previous wavefunction can be recast in the form 



i<jGA 



i<jeB 



(107) 



where we made explicit the Jastrow factor in each cluster and A refers to the antisymmetrization over all the possible clusterings 
in the sets A and B. In the spirit of the quantum Hall states showing a clustering into two sets - the main example being the 
Moore and Read (MR) Pfaffian state |60, 61] - i^c is characterized by the presence of quasi-hole excitations corresponding to 
half a quantum flux (and effective charge 1, since = 2). These excitations must appear in pairs: it is however helpful to analyze 
the two possible wavefunctions they can assume. Such wavefunctions are related to quasi-holes in the two different clusters, 
given by 



<Ji (C) = A 
02 (C) = A 



n ^i.'^ n - n -on ^2,/ n - ^^■) 

A A A B B 

n si,k n - n ^2./ n - ^^■) n - 



B B 



2N 

e 1;...;) 



2Af 

e 1;...;) 



(108) 
(109) 



These quasi-holes obey a fermionic statistics: once two of them of the same kind are exchanged, the wavefunction acquires a tt 
pha se. However, it is interesting to notice that they show the same fusion rules of the Ising model with defined fermionic parity 
06211 (characterizing the MR state i60ll6lll ) once one defines a third bosonic excitation given by the fusion ij: = ai x a2- From 
the physical point of view, if it is possible to obtain linear superpositions of ci and a2 through the interplay between repulsive 
potentials for the |t) and \l) components, then the so-obtained (non-Abelian) excitations could present interesting features, a 
point which certainly deserves further investigations. 

In fermionic systems at the degeneracy point = 3B, the ground-state ilc, characterized by the filling factor = 2, is the 
highest density state obtained with atoms in the Hilbert space spanned by the deformed Landau level Xi and X2 ■ Moreover, ilc 
minimizes the term in the Hamiltonian proportional to the angular momentum. Nevertheless, in the study of rotating ultracold 
atomic gases, it is interesting to analyze what happens varying the filling factor, since, in general, such systems present non-trivial 
phase diagrams as a function of IM HI . As we have already shown, the double degeneracy at this particular value of q 
provides in a natural way a clustering of the atoms into two sets in order to minimize \J\. Each atom can assume a wavefunction 
which is a superposition of states in Xi and in X2 and transitions from one to the other are possible: therefore, also for smaller 
values of the filling factor, we are driven to consider deformed ground-states showing a pairing among tha atoms that are similar 
to the ones usually considered in the study of fractional quantum Hall effect |64]. 

Let us consider first the filling factor = 1: in this case a paired state is built by favouring the creation of coupled atoms in 
the antisymmetric state obtained by applying the operator {Gi,iG2,j — G2,iGi,j) to the pair of atoms (i, j) in the limit Zi — > zj. 
The corresponding wavefunction for 2N particles reads 



ilm = Hf (Gi,iG2,j — G2,iGi, 



1 



2N 



Y[{zi - Zj) 



IN 

e ^ - lU ... i) 



(110) 



where Hf indicates the Haffnian, which is a symmetric version of the Pfaffian defined for a symmetric matrix My = M^^. 



Hf {M) = 



'(4)---^ct(2A'-1),(t(2JV) 



(111) 



(the sum is over all the permutation of the indices). 

The wavefunction finf is antisymmetric over all the atoms because it is composed by the symmetric Haffnian and by the 
antisymmetric Jastrow factor. This implies that intra-species interactions give a zero contribution to the energy and finf can 
be considered a ground-state since each atom lies in a superposition of states of Xi and X2 ■ However this wavefunction is 
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not vanishing for Zi Zj because of the components with different spin, therefore Jlnf is not in general a ground-state for 
inter-species repulsive interactions that require higher powers of the Jastrow factor to give a null contribution. 

For 1/ = 1/2 there are two possible antisymmetric paired states that have been widely analyzed in the literature. The first 
one corresponds to a deformed MR Pfaffian state |60, 61] and the second corresponds to a deformed Haldane-Rezayi state 
1160, 65J. The deformed MR Pfaffian state can be described by an effective p-wave pairing [64J obtained by applying the operator 
{Gi,iG2,j + G2,iGi,j) which favours a symmetric state (with respect to Xi X2) for the pair when Zi — > Zj. The 
corresponding wavefunction is ll44ll 

- Pf {{Gl,^G2,o + G2^£l,o) H " ^j)' ^^^^^ - ^) ^l^^) 

Zi Zr 



where Pf is the Pfaffian operator. This wavefunction is antisymmetric, therefore intra-species interactions give a null contribu- 
tion and VLmr can be considered a ground-state since each atom lies in a superposition of states in Xi and X2 ■ This state shares 
all the main characteristics of the Moore and Read wavefunctions |60,], and, in particular, its excitations are non-Abelian Ising 
anyons, as shown in ll6lll where an analogous wavefunction is analyzed. VImr can be mapped into the usual spinless MR state 
'^MR ll60l I61I1 in a way which is similar to equation (|95] l for the ground-state outside the degeneracy points. Since, for every 
factor in the Pfaffian, one atom is in a state in Xi and the other in X2^ the norm of 0,mr is obtained from the one of the Pfaffian 
state MR just by multiplying it by a constant value for each pair of atoms: 

{^mr\ ^mr) = (|c^4|'+8S|c^^i|')'^(8B|ct,2|'+2(8B)2|c^^2|')'^(*Mfl| ^a/t?,) • (113) 



This constant value can be calculated in a way similar to equation ( |96] l and it is an effect of the clustering characterizing the state 
^Imr- Equation ( II 13b guarantees that also the statistics of the excitations (a and (i, of the kind 

2N 

n t/- /■\ T3f fSi,iG2,iizt~Ca){zj -Cb)+i^ j\Yj , s2 -f 1:1^.1' , 

i^Mfl (Ca,C&) = Pf ^ j[[(Zi-Zj) e • \U ■■■ i) (114) 

can be described by Ising anyons as in the case of the Moore and Read state, since the Berry and monodromy phases acquired 
in the exchange of two excitations coincide. 

The other paired ground-state at ly — 1/2 is a deformed Haldane-Rezayi state that can be obtained through the introduction 
of the antisymmetric operator {Gi,iG2,j — G2,iGi.j)- 

( 1 \ 

^HR=Vi\(9lA.J-G2,^Gl,i)- 72 n(-^»--^j)'e ' * ^ (115) 

This state has a total angular momentum | J| which is lower than the one of ^mr. so that, in principle, it is energetically favoured 
if we consider the single-particle Hamiltonian ( [36] l. However the Haldane-Rezayi state represents the critical point between 
a weak and a strong coupling phase in a fermionic system with an effective c?-wave pairing |[64t], therefore it is considered 
to be a gapless state whose excitations are described by a non-unitary conformal field theory \W\ which is unfit to define an 
incompressible state. Moreover we expect that JIht? is more influenced by the presence of a weak inter-species interaction than 

^MR- 

To conclude the discussion about the (free) fermionic gases it is worth noticing that all the states presented can be retrieved 
also for the generic degeneracy point between states in x^ and Xn+i substituting the operators G\ and G2 with Gn and Gn+i- 
This is true in the case of free fermions (i.e., with go — 0), whereas the introduction of strong interactions, such as an inter- 
species contact repulsion, brings to different scenarios presenting deformed Halperin states, similar to the one we will present 
for the bosonic gases, whose exponents and filling factors depend also on the value of n. 



B. Bosonic gases 

The analysis of the degenerate points can be applied also to bosonic gases. As described in Section IV A| the state 
(llOlb is a ground-state of both the intra-species and the inter-species contact interactions for every m > 2n, therefore 
is the bosonic wavefunction with zero interaction energy minimizing the polynomial order In the point = 3B one has a 
superposition of states in Xi and in X2 and, as in the previous case, we can describe the related multiparticle ground-state 
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through the clustering in two corresponding subsets A and B. The resulting ground-state is 



n^i ^ n 

feeA i<jeA 



n^- n 

;es i<jeB 



keA.ieB 



e • \ii ... i) 



S 



keA 



leB 



KjeB 



A 



(4) 
2JV 



(116) 



where the symmetrization over all the atoms is necessary to have a bosonic state and A*^"*^ is the generaUzed Laughlin state ( |94] l. 
The exponents of the Jastrow factors are defined considering that the lowest polynomial order term for an atom in A is given 
by the first derivative included in Qi, whereas in B by the second derivative in Q2. Therefore fif, vanishes at least as (z^ — Zj) 
whenever Zi — > Zj for each pair of atoms, and the intra-species interaction energy is zero, fif, is characterized by a filling factor 
1/ = 2/9 and it is built by the symmetrization of the Halperin state (4, 6, 4). The interaction term between atoms in different 
clusters is determined in order to satisfy the zero interaction energy constraint 



keA,ieB 



The ilfc filling factor, = 2/9, is therefore the highest possible filling factor that guarantees a null interaction energy for bosons 
at the point = 3B. 

We conclude this Section observing that considering different degeneracy points, = (2n + 1)B, the ground states for intra- 
species repulsions are defined as deformed Halperin states (2n + 2,2n + 4, 2n + 2) characterized by lower filUng factors. 



C. Crossover at a degeneracy point 



So far we considered the non-Abelian component of the potential at the exact value q^ = (2n + l)B which is characterized by 
a perfect degeneracy of the deformed Landau levels Xn ^i^^ Xn+i- Then the population of atoms must be equally distributed into 
these levels in order to minimize the total angular momentum of the system and thus its total energy. However, if the parameter 
q^ is slightly detuned from these degeneracy points, the deformed Landau level with a lower energy will present an increase in 
population balancing the small energy difference arising between the two state famiUes. 

Let us consider, in particular, a system showing a filling factor = 1 /m„ on the n}^ Landau level, and let us suppose that 

To fill the gap between Xn ^i^^ Xn+i there must be an imbalance 



is sHghtly higher than {2n + l)B, so that > £„+i.,„. 



M = N, 



Nn between the population of atoms in Xn+i Xn ' ^^^^ that 



Considering the energy eigenvalues given by Equation ( 11021 ), one obtains 



Mto„+iA sa 5q^ 



■ e. 



ri+i,o; 



dg2 



q^ = (2n+l)B 



1 



1 



4n + 3 An+l 



(117) 



where Sq^ is the displacement of q^ from the degeneracy point: we assumed A <C i? in the second approximation done in 
Equation ( 1117b . In order to obtain a system described by the deformed Hall states defined above for the degeneracy point at 
q^ — 3B, the imbalance M must be negligible with respect to the total number of atoms N and, in particular, one obtains the 
following degeneracy condition on the displacement Sq^ around 3B: 



12 I. 21 

— \6q^\ « -A < 
35 ' ' u 



9 



73 B w 0.46B, 



(118) 



where v is the total filling factor of the system and the further constraint for NA is derived from the energy difference with the 
third Landau level, X3 , in such a way that all the atoms lie only in the two lower Landau levels. If Sq^ exceeds the limit ( II 181 ). 
the imbalance M grows up until only the population in the lower Landau level is left. Therefore the displacement 5q^ drives a 
crossover between different regimes characterized by wavefunctions in Xn ^"d Xn+i such as the deformed Laughlin states in 
(llOlb and only in the regime defined by the condition dl 18b Hall states presenting particles in both the deformed Landau levels 
can be present. 
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Vm. CONCLUSIONS 

We studied two-component ultracold bosonic and fermionic atomic gases in U{2) non-Abelian potentials. We focused our 
attention to gauge potentials having, in the Abelian limit, doubly degenerate Landau levels deriving from an Abelian magnetic 
field equal for both the components. We investigated the effect of general homogeneous non-Abelian terms, discussing the 
conditions under which the structure of degenerate Landau levels is preserved (even if the spin degeneracy is split by the 
coupling between the internal states realized by the non-Abelian term). We argued that there are only three classes of (quadratic) 
Hamiltonians preserving the degeneracy of the Landau levels and giving rise to an analytically defined Landau level structure, 
where the eigenstates are expressed as a finite linear combination of the eigenstates of a particle in a magnetic field. The first 
class corresponds to an Abelian U{1) x U{1) gauge potential and it refers to uncoupled internal states, whereas the other two 
are characterized by a truly non-Abelian U{2) gauge potential and correspond to different kinds of Jaynes-Cummings models. 

Focusing on one of these Jaynes-Cummings classes which is gauge equivalent both to a Rashba and to a DresseUiaus spin- 
orbit coupling, we determined the single-particle energy spectrum and we discussed the parameters of the Rabi pulses needed for 
the physical implementation of such gauge potentials. It is found that the single-particle energy levels assume the same values 
in a series of degeneracy points. The corresponding deformed quantum Hall states for fermions and bosons (in the presence 
of strong intra-species interaction) have been determined through a mapping to the usual quantum Hall Laughlin states. Such 
mapping preserves the Berry phases characterizing the exchanges of quasi-holes and, therefore, the Abelian excitation statistics. 
Far from the degeneracy points, deformed Laughlin states are found; whereas, at the crossing points of the lowest Landau levels, 
ground-states with non-Abelian excitations emerge at different fillings. A detailed discussion of the properties of the resulting 
deformed Moore-Read states and of the state ilc (minimizing the total angular momentum) was presented; besides, the crossover 
arising at the degeneracy points was investigated, estimating the range of the Hamiltonian parameters for the validity of these 
many-body ground states. 

The effect of an angular moment term on the stability of the Landau level structure and of the deformed Laughlin states was 
analyzed and we provided an estimate of the values of the angular momentum such that the Landau level description of the 
multiparticle states remains accurate. We also gave results for U (3) gauge potentials acting on a three-component gas, pointing 
out that it is possible to have lines of degeneracy points of the Landau levels and triple degeneracies in the spectrum. 

It would be interesting in the future to study the transport properties of the quantum Hall states found for ultracold atoms in 
U{2) artificial gauge potentials preserving the Landau levels, especially in view of the experimental study of the signatures of 
the different Hall states. From this point of view, we mention that it could be useful to study setups in which the Abelian terms 
are different for the two components. Another related issue, important for the experimental detection of non-Abelian states, 
would be also the investigation of the physical addressability of such states, e.g. how the filling of these states depends on the 
experimental parameters. Equally important would be the estimation of Haldane pseudopotentials relevant to address the effect 
of more general interacting terms like inter-species repulsions and dipolar long-range interactions. 
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